Probability Estimating Apparatus and Metliod for Peak-to-Peak Clock Skews 



This application claims the benefit of U.S. Provisional Application No. 
60/277,251 filed March 20, 2001. 

Field of the Invention 

This invention relates to a probability estimating apparatus and method for peak- 
to-peak clock skews for testing the clock skews among a plurality of clock signals 
distributed by a clock distributing circuit, and for estimating the generation probability of 
the peak-to-peak value or peak value of the clock skews. 

Background of the Invention 

Traditionally, the clock skews, as shown in Figure 1, are statistically estimated by 
using a time interval analyzer or a fi-equency counter. The time interval analyzer 
measures the time difference of zero-crossing points between the clock signals under test 
and a reference clock signal, then measures the fluctuations of the time difference 
through histogram analysis. An example of such a clock skew measurement using a time 
interval analyzer is described in Jitter Analysis Clock Solutions, by Wavecrest 
Corporation, 1998. 

When evaluating reliabiUty of a microprocessor, for example, it is effective to 
detennme whether the probability of peak-to-peak clock skews in the clocks distributed 
within the microprocessor exceeds a predetermined value, or to confirm that the 
generation probability of the peak-to-peak value in the clock skews will not exceed the 
predetermined value. However, since the generation probability of the peak-to-peak 
value of the clock skews have never been theoretically analyzed, the traditional method 
requires an enormous amount of data to estimate the generation probability of the peak- 
to-peak value of the clock skews, and thus, requiring a large amount of time. 

Summary of the Invention 

It is, therefore, an object of the present invention to provide a probability 
estimating apparatus and method for peak-to-peak clock skews for estimating the 
probability of the peak-to-peak value among the clock signals under test to see whether it 
exceeds the predetermined value, as well as estimating the generation probability of the 
peak-to-peak value in the clock skews in a much shorter time. 
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The above object is achieved by the probability estimating apparatus for peak^^ 
peak clock skews of the present invention which is characterized by having a clock skew 
estimator for estimating the clock skew sequences between a plurality of clock signals 
under test, and a probability estimator for determining and outputting the generation 
5 probability of the peak-to-peak value of the clock skews mnong the clock signals under 
test based on the clock skew sequences. 

Further, in the probability estimating apparatus for peak-to-peak clock skews, 
preferably, the probability estimator determines and outputs the generation probability of 
llie peak-to-peak value of the clock skews among the clock signals under test based on 
10 the clock skew sequences. 
CJI Further, in the probability estimating apparatus for peak-to-peak clock skews, 

preferably, the probability estimator is comprised of an RMS (root mean square) detector 
for determining the RMS value of the clock skew sequence data suppUed thereto, a 
|1 memory for storing a predetermined value, and a probability calculator for determining 

15 attid outputting the probability of the peak-to-peak clock skews among the signals under 
C«i test to determine whether the probability exceeds the predetermined value based on the 

pj predetermined value and llie RMS value. 

II Further, in flie probability estimating apparatus for peak-to-peak clock skews, 

preferably, the probability estimator is comprised of an RMS detector for determining the 
20 RMS value of the clock skew sequence data supplied thereto, a peak-to-peak detector for 
calculating the difference between the maximum and minimum values based on the clock 
skew sequence data to determine the peak-to-peak value, and a probability calculator for 
obtaining and outputting the probability of the clock skews among the clock signals 
under test to determine whether the probability exceeds the peak-to-peak value based on 
25 the peak-to-peak value and the RMS value of those clock skew sequence data. 

Further, in the probability estimating apparatus for peak-to-peak clock skews, 
preferably, the clock skew estimator is furnished with a timing jitter estimator for 
estimating the timing jitter sequences of a plurality of clock signals under test, and a 
clock skew calculator for receiving a plurality of those timing jitter sequences as inputs 
30 and calculating the timing difference sequences among the timing jitter sequences to 
output the clock skew sequences. 
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Further, in the probabiUty estimating apparatus for peak-to-peak clock skews, 
preferably, the clock skew estimator is furnished with a second clock skew calculator for 
receiving the clock skew sequences as inputs to determine the difference among the 
plurality of those clock skew sequences. 
5 Further, in the probability estimating apparatus for peak-to-peak clock skews, 

preferably, the clock skew estimator is furnished with a frequency multiplier for 
receiving the timing jitter sequences as inputs to produce timing jitter sequences by 
multiplying the frequency of the clock signals under test. 

Further, in the probabiUty estimating apparatus for peak-to-peak clock skews, 
10 preferably, the clock skew estimator is furnished witii a deterministic clock skew 

estimator for estimating the timing errors among ideal clock edges of a plurality of clock 
signals imder test to output the deterministic components of those clock skews. 

Further, in the probability estimating apparatus for peak-to-peak clock skews, 
preferably, the timing jitter estimator is comprised of an analytic signal transformer for 
U^. -. 15 transforaiing the clock signals under test into analytic signals of complex nimiber, an 
instantaneous phase estimator for determining an instantaneous phase of the analytic 
ll^j signals, a linear phase remover for removing the linear phase from the instantaneous 

I J phase to obtain instantaneous phase noise, and a zero-crossing resampler for receiving the 

instantaneous phase noise as inputs and resampling only the instantaneous phase noise 
20 data closest to the zero-crossing timings of the real number in the analytic signals to 
output the timing jitter sequences. 

Further, in the probability estimating apparatus for peak-to-peak clock skews, 
preferably, the analytic signal transformer is comprised of a band-pass filter for receiving 
the clock signals under test and extracting only the components closest to the 
25 fimdamental frequency from the clock signals under test to band-limit the clock signals 
under test, and a Hilbert transformer for Hilbert-transforming the output signals of the 
band-pass filter to generate a Hilbert converted pair of the clock signals. 

Further, in the probabiUty estimating apparatus for peak-to-peak clock skews, 
preferably, the analytic signal transformer is comprised of a time domain to frequency 
30 domain transformer for receiving the clock signals under test and transforming the clock 
signals into both-side spectra signals in a frequency domain, a bandwidth limits for 
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extracting only the components closest to the positive fundamental frequency from the 
both-side spectra signals, and a frequency domain to time domain transformer for 
transforming the outputs of the bandwidth limiter back to the time domain signals. 
Further, in the probability estimating apparatus for peak-to-peak clock skews, 
5 preferably, the analytic signal transformer is comprised of a buffer memory for storing 
the clock signals under test, means for sequentially extracting the signals from the buffer 
memory while overlapping a part of the extract signals with the ones previously 
extracted, means for multiplying a window ftmction by each extracted part of the signal, 
means for transforming each multipUed part of the signal into the both-side spectra 
10 signals in the frequency domain, a bandwidth limiter for extracting only the components 
closest to the positive fundamental frequency of the clock signals under test from the 
both-side spectra signals transformed into the frequency domain, meaiis for transform 
M the outputs of the band-pass filter back to time domain signals, and means for multiplying 



an inverse of the window function by the signals transformed to the time domain to 
15 obtain band-limited analytic signals. 

Further, in the probability estimating apparatus for peak-to-peak clock skews, 
preferably, the clock skew estimator is furnished with an AD converter for receiving the 
clock signals under test as inputs and digitizing the analog signals to covert to digital 
LJ signals. 

20 Further, in the probability estimating apparatus for peak-to-peak clock skews, 

preferably, the clock skew estimator is furnished with a waveform clipper for receiving 
the clock signals under test as inputs and removing the amplitude modulation components 
to extract only the phase modulation components of tiie clock signals under test 

Further, in the probability estimating appdmtas for peak-to-peak clock skews, 

25 preferably, the analytic signal transformer has an adjustable pass band for the clock 
signals under test. 

Further, in the probability estimating apparatus for peak-to-peak clock skews, 
preferably, the timing jitter estimator is furnished with a low frequency component 
remover for receiving the instantaneous phase noise as inputs and removing the low 
30 frequency components from the instantaneous phase noise before providing the 
instantaneous phase nose to the zero-crossing resampler. 
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Finther, the above object is achieved by the probabihty estimating method for 
peak-to-peak clock skews of the present invention characterized by having a step of 
estimating the clock skew sequences among a plurality of clock signals under test, and a 
step of determining and outputting the generation probability of the peak-to-peak value of 
5 the clock skews among the clock signals under test based on the clock skew sequences. 

Further, the object is achieved by the probability estimating method for peak-to- 
peak clock skews characterized by having a step of estimating the clock skew sequences 
among a plurality of clock signals under test, and a step of determining and outputting the 
generation probability of the peak value of the clock skews among the clock signals 
10 under test based on the clock skew sequences. 

Further, in the above probability estimating method for peak-to-peak clock skews, 

in ■ 

5k the step for determining the generation probability of the clock skew peak value is 

•%i . . 

W preferably comprised of a step of determining an RMS (root mean square) value of the 

p clock skew sequence data supplied thereto, and a step of determining the probability for 

15 the peak-to-peak clock skews among the clock signals under test which exceeds the 

p predetermined value based on tiie predetermined value and the RMS value of the clock 

W 

1 1 skew sequence data. 

|;;! Fxirther, in tiie above probability estimating method for peak-to-peak clock skews, 

|,| the step for determining the generation probabihty of the clock skew peak value is 

20 preferably comprised of a step of determining the RMS value of the clock skew sequence 
data supplied thereto, a step of calculating the difference between the maximum and 
minimum values of the clock skew sequence data to determine the peak-to-peak value, 
and a step of determining the probabihty of the clock skews among the clock signals 
imder test that exceeds the peak-to-peak value based on the peak-to-peak value and the 
25 RMS value of the clock skew sequence data. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
step for estimating the clock skew sequence data is preferably comprised of a step of 
estimating the timing jitter sequences of a plurality of clock signals under test, and a step 
of receiving the plurahty of timing jitter sequences and calculating differences among the 
30 timing jitter sequences to estimate the clock skew sequences. 
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Further, in the probabiUty estimating method for peak-to-peak clock skews, the 
step for estimating the clock skew sequences is preferably comprised of a step of 
receiving the above clock skew sequences and determining the differences among the 
plurality of the clock skew sequences to estimate the clock skew sequences. 
5 Further, in the probability estimating method for peak-to-peak clock skews, flie 

step for estimating the clock skew sequences is preferably comprised of a step of 
receiving the timing jitter sequences and estimating the timing jitter sequences by 
multiplying the frequency of the clock signals under test. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
10 step for estimating the clock skew sequences is preferably comprised of a step of 

estimating the timing errors among the ideal clock edges of a plurality of clock signals 
W under test to estimate the deterministic components of the clock skews, 

w Further, in the probability estimating method for peak-to-peak clock skews, the 

step for estimating the timing jitter sequences is preferably comprised of a step of 

W ;.. 

II 15 transforming the clock signals under test into analytic signals of complex numbers, a step 

m: 

of determining the instantaneous phase of the clock signals imder test based on the 
analytic signals, a step of removing the linear phase from the instantaneous phase to 
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f !| estimate the instantaneous phase noise, and a step of receiAdng the instantaneous phase 

^ ^ noise and resampling only the instantaneous phase noise data closest to the zero-crossing 

20 timings of real number in the analytic signals to output the timing jitter sequences. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
step for transforming the signals under test into anal5^c signals is preferably comprised 
of a step of extracting only the components closest to the fundamental frequency of the 
clock signals xmder test to band-limit the clock signals under test, and a step of Hilbert- 
25 transforming the output signals from the band-pass filter to generate a Hilbert 
transformed pair of the input signals. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
step for transforming the clock signals under test into analytic signals is preferably 
comprised of a step of transforming the clock signals under test into both-side spectra 
30 signals in a frequency domain, a step of extracting only the components closest to the 
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positive fundamental frequency from the both-side spectra signals, and a step of 
transfomiing the outputs of the bandwidth limiter back to time domain signals. 

Further, in the probabihty estimating method for peak4o-peak clock skews, 
step for transforming the clock signals under test into analytic signals is preferably 
5 comprised of a step of collecting the clock signals under test in a buffer memory, a step 
of sequentially extracting the signals from the buffer memory while overlapping a part of 
the extracted signals with the ones previously extracted, a step of multiplying a window 
fttaction by each extracted part of the signal, a step of transforming each multiplied part 
of the signal into both-side spectra signals in the frequency domain, a step of extracting 
10 only the components closest to the positive ftindamental frequency of the clock signals 
under test from the both-side spectra signals transformed into tihe frequency domain, a 

CI 

g-l step of transforming the band-limited spectra signals into time domain signals, and a step 



of multiplying an inverse window ftinction by the signals transformed into the time 
Ml domain to obtain band-limited analytic signals. 

^ 15 Further, in the probability estimating method for peak-to-peak clock skews, the 

St step for estimating the deterministic components of the clock skews among the clock 

signals under test is preferably comprised of a step of determining the deterministic 



components of the clock skews by receiving the linear instantaneous phase of a plurality 
of clock signals under test and determining the difference between the initial phase angles 

20 of the linear instantaneous phase. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
step for estimating the deterministic components of the clock skews among the clock 
signals under test is preferably comprised of a step of estimating the clock edges 
corresponding to the clock signals and determining the offset value of the clock edges by 

25 receiving the timing jitter sequences of a plurality of clock signals under test and 
determining the correlation among the timing jitter sequences. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
step for estimating the deterministic components of the clock skews among the clock 
signals under test is preferably comprised of a step of determining the deterministic 

30 components of the clock skews by receiving the plurality of clock signals under test and 
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determining the average of the zero-crossing timing errors among the clock signals under 
test. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
step for estimating the clock skew sequences is preferably comprised of a step of 
conducting waveform clipping on the clock signals under test to remove the amplitude 
modulation components, thereby extracting only the phase modulation components from 
the clock signals under test. 

Further, in the probability estimating method for peak-to-peak clock skews, the 
step for estimating the timing jitters is preferably comprised of a step of receiving the 
instantaneous phase noise and removing the low frequency components from the 
instantaneous phase noise. 

Brief Description of the Drawings 

Figure 1 is a diagram showing an example of clock skew measurement using a 
time interval analyzer. 

Figure 2A is a schematic diagram showing a clock network for distributing clock 
signals, and Figure 2B is a schematic diagram showing the timing relationships among 
the clock skews. 

Figure 3 is a schematic diagram showing the timing relationships between the 
timing jitters and the clock skews. 

Figure 4 is a diagram showing an example of waveforms of the clock signal under 

test. 

Figure 5A is a diagram showing an example of timing jitter A^[n] of the clock 
signal Xj(t) under test, and Figure 5B is a diagram showing an example of timing jitter 
A<|)^[n] of the clock signal Xk(t) under test. 

Figure 6 is a diagram showing an example of clock skew between the clock 
signals imder test measured by the clock skew measurement method of the present 
invention. 

Figure 7 is a diagram showing an example of histogram of the clock skew 
between the clock signals under test measured by the clock skew measurement method of 
the present invention. 
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Figure 8 is a schematic diagram showing a clock distribution network with 
different clock domains. 

Figure 9 is a schematic diagram showing a principle of clock skew measurement 
of the present invention involving frequency multiplication. 

Figure 10 is a diagram showing an example of power spectra of a clock signal in a 
microprocessor. 

Figure 11 is a diagram showing an example of timing jitter histogram of a clock 
signal in a microprocessor. 

Figure 12 is a diagram showing a probability density fimction of the Rayleigh 
distribution. 

Figure 13 is a diagram showing a right-tail probability of the Rayleigh 
distribution. 

Figure 14 is a diagram showing a relationship between the number of zero- 
crossing samples and the observed peak-to-peak clock skew values. 

Figure 15 is a diagram showing an example of waveform of the clock signal under 

test. 

Figure 16 is a diagram showing an example of waveform of an analytic signal of 
the clock signal under test. 

Figure 17 is a diagram showing an example of waveform of the instantaneous 
phase of the clock signal under test. 

Figure 18 is a diagram showing an example of waveform of the instantaneous 
phase noise of the clock signal under test. 

Figure 1 9 is a diagram showing an example of waveform of the timing jitter of the 
clock signal under test. 

Figure 20 is a diagram showing an example of waveform of the of clock signal 
under test. 

Figure 2 1 is a diagram showing an example of waveform of the Hilbert 
transformed analytic signal. 

Figure 22 is a diagram showing an example of instantaneous phase signal having 
a discontinuous instantaneous phase. 



Figure 23 is a diagram showing an example of continuous instantaneous phase 
signal that has been unwrapped. 

Figure 24 is a diagram showing an example of digitized clock signal under test. 

Figure 25 is a diagram showing an example of both-side power spectra of the 
clock signal under test obtained through FFT. 

Figure 26 is a diagram showing an example of one side of the band-limited power 
spectra of the clock signal. 

Figure 27 is a diagram showing an example of band limited analytic signal 
derived through inverse FFT. 

Figure 28 is a diagram showing an example of adaptive zero-crossing points 
approximation of the clock signal under test. 

Figure 29 is a diagram showing an example of clock signal under test having AM 
components. 

Figure 30 is a diagram showing an example of clock signal under test without 
having AM components. 

Figure 3 1 is a diagram showing an example of structure in the ^paratus for 
estunating the probability of peak-to-peak clock skews in the present invention. 

Figure 32 is a flow chart showing an example of method for estimating the 
probabihty of peak-to-peak clock skews in the present invention. 

Figure 33 is a diagram showing another example of structure in the apparatus for 
estimating the probability of peak-to-peak clock skews in the present invention. 

Figure 34 is a flow chart showing another example of method for estimating the 
probability of peak-to-peak clock skews in the present invention. 

Figure 35 is a diagram showing an example of structure in the clock skew 
estimator in the present invention. 

Figure 36 is a flow chart showing an example of clock skew estimation method in 
the present invention. 

Figure 37 is a diagram showing an example of structure in the deterministic clock 
skew estimator in the present invention. 

Figure 38 is a flow chart showing an example of deterministic clock skew 
estimation method in the present invention. 
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Figure 39 is a diagram showing another example of structure in the clock skew 
estimator in the present invention. 

Figure 40 is a flow chart showing another example of clock skew estimation 
method in the present invention. 

Figure 41 is a diagram showing a further example of structure in the clock skew 
estimator in the present invention. 

Figure 42 is a flow chart showing a further example of clock skew estimation 
method in the present invention. 

Figure 43 is a diagram showing an example of structure in the timing jitter 
estunator used in the clock skew estimator in the present invention. 

Figure 44 is a flow chart showing an example of timing jitter estimating method 
used in the clock skew estimating method in the present invention. 

Figure 45 is a diagram showing an example of structure in the analytic signal 
transformer used in the clock skew estimator in the present invention. 

Figure 46 is a flow chart showing an example of analytic signal transformation 
method used in the clock skew estimating method in the present invention. 

Figure 47 is a diagram showing another example of structure in the analj^c signal 
transformer used in the clock skew estimator in the present invention. 

Figure 48 is a flow chart showing another example of analytic signal 
transformation method used in the clock skew estimating method in the present 
invention. 

Figure 49 is a diagram showing a further example of structure in the analytic 
signal transformer used in the clock skew estimator in the present invention. 

Figure 50 is a flow ch^ showing a further example of analytic signal 
transformation method used in the clock skew estimating method in the present 
invention. 

Figure 51 is a diagram showing a further example of structure in the clock skew 
estimator in the present invention. 

Figure 52 is a flow chart showing a further example of clock skew estimation 
method in the present invention. 



11 



Figure 53 is a diagram showing a further example of structure in the clock skew 
estimator in the present invention. 

Figure 54 is a flow chart showing a further example of clock skew estimation 
method in the present invention. 

Figure 55 is a diagram showing a further example of structure in the timing jitter 
estimator used in the clock skew estimator in the present invention. 

Figure 56 is a flow chart showing a ftirther example of timing jitter estimating 
method used in the clock skew estimating method in the present invention. 

Detailed Description of the Invention 

The structure and operations of the present invention will be described in detail 
below. In the following, MPU clock signals are used as signals to be tested for the 
purpose of explanation of the invention. 

First, the clock skew will be defined. As shown in Figure 2A, where, for 
example, the clock source of the clock distribution network is used as a reference point, 
the clock skew is defined as a difference between delay times T^cd and T^cd before clock 
signals CLKj and CLE^c arrive at registers R, and R^c. 

TiLinT)=TUnT)-TUnT) 

(1) 

Figure 2B shows the timings of the clock skews. Here, T is a fundamental clock period 
of the clock signals under test. 

The rising edge times of the clock signal CLKg, CLK^, and CLECk are denoted by 
t®cd, f ed, and ^cd respectively. Also, when the ideal clock edge time (the jitter-free clock 
edge time) of each clock signal CLKg, CLKj, and CLKk is denoted by (nT)g, (nT)j, and 
(nT)k, the delay times Vcd and T^'cd will be respectively expressed as follows: 

= k(«r)-(nr)J-b(„r)-(«r)J+^„r), -(„r)J 




[sec] (2) 
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~ ^ Skew 



In J 



[sec] (3) 

The following eqxiations express the respective tune differences between the ideal clock 
edge times of CLKj, CLKk and the ideal clock edge time CLKg, and correspond to the 
deterministic components (deterministic clock skew values) of the clock skew determined 
by signal paths. 



[sec] 



(4) 



[sec] (5) 

In addition, A<|>8[n](Tg/27t) (= t^cd (nT) - (nT)g), M^[n]{yin) (= iJcd (nT) - (nT)j), and 
A(|)''[n](Tk/27t) (= t^cd (nT) - (nT)k) express the timing jittw sequences (unit in second) of 
the clock signals CLKg, CLKj, and CLKk respectively. When equations (2) and (3) are 
substituted into equation (1), the clock skew T''''skew in CLKj and CLKk will be estimated 
as follows: 







\-l^'[n] 








-t4'[n] 





^ Skew ^ 



(T \ It 
2nJ ^ ^ \2n 



[sec] (6) 

the second member in the equation (6) as noted below corresponds to the random 
variation (random component) of the clock skew based on the timing jitters in each clock 
signal. 
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In other words, this clock skew estimation method can determine the random distribution 
of a clock skew by combining the amount of difference between the clock edge of each 
clock signal and the ideal clock edge, namely, the timing jitters of each clock signal. 
Here, the fundamental periods of the generally distributed clock signals CLKj and CLKk 
are identical to one another (Tj = Tk). Figure 3 shows the relationship between the timing 
jitter and clock skew. 

The following equation expresses the difference between the rising edge times of 
the ideal clocks of the clock CLKj and CLKk, and are the deterministic components of the 
clock skew determined by the signal paths of the clock distribution network: 

[sec] (7) 

The deterministic clock skew value x^'^skew, for example, can be determined by the 
instantaneous phases of two signals CLKj, CLKk under test and the difference between 
their linear phase components in the instantaneous phases. The fundamental cosine wave 
components of the signals CLKj and CLKk are as follow: 



Xj{t)=^Aj cos{i^' (t))=Aj cos ^t+^i -A(|)^(0 



xM^A^ cosi^'{t))=A^ cos ^r+<|)o -^'(t) 



. T 



(8) 



(9) 



Here, the instantaneous phases of Xj(t) and Xk(t) are respectively expressed by the sum of 
the instantaneous linear phase component 2nt/TL having the fimdamental period Tl (L = 
j, k), initial phase angle ^\ (L = j, k), and instantaneous phase noise component A^\t) (L 



j 



[rad] 



(10) 



[rad] 



(11) 



14 



The estimation method for the instantaneous phase in the clock signals will be explained 
later When A^(t) = 0 is used in equations (10) and (1 1), the instantaneous linear phase in 
the jitter-free clock signals are obtained as follows: 



linear 



[rad] 



(12) 



^linear 



[rad] (13) 

Here, the ideal rising edge times t = (nT)j, and t = (nT)k of the signals CLKj and CLKk 
are the times where the instantaneous linear phase becomes (Irm - n/2), and have the 
following relationsiiip from equations (12) and (13): 



2n 



[sec] 



[sec] 



(14) 



(15) 



Therefore, per equation (7), the following deterministic clock skew value is obtained: 



n 



2nn <|)o 



1^ 

27t 



271 



= <l>0^ 



2n 



-<!>o 



271 



t4/-*.*)|- 



[sec] (16) 

In general, the fundamental periods of the distributed clock signals CLKj and CLKk are 
the same with one another (T, = Tk). In other words, the deterministic clock skew value 
between the two signals under test can be determined as a difference between the initial 
phase angles in the instantaneous linear phase of those two signals under test. 

Here, the initial phase angle (|)o of the signals under test can be obtained by 
conducting a linear line fitting for the instantaneous phase waveform data (|)(k) based on 



15 



the minimum square method and choosing the formula below in such a way that 
becomes the minimum: 



Here, the initial phase angle to be determined is expressed as follows: 
2N(2N+l)f^4k)-6f^k4k) 



(17) 



(18) 



The estimation of parameters based on the linear line fitting is described, for example, 
page 362 of J. S. Bendat and A. G. Piersoi "Random Data: Analysis and Measurement 
Procedure'', 2""^ edition, published by John Wiley & Sons, Inc. in 1986. 

Further, the initial phase angle of signal x(t) under test can be obtained by 
conducting a cosine wave fitting for the clock waveform data x(k) or fundamental sine 
wave component based on the minimum square method in the following equation in such 
a way that becomes the minimum through the maximum likelihood estimation method: 



2 x{kyAcos^k + ^, 



Here, the initial phase angle to be determined is expressed as follows: 



(19) 



^0-- arctan 



k=l i 

^x{k)cos~k 



(20) 



The estimation of the parameter based on the maximum likelihood estimation is 
described, for example, pages 167-172 of S. M. Kay ''Fundamentals of Statistical Signal 
Processing: Estimation Theory"\ pubUshed by Prentice-Hall, Inc. in 1993. 

In the foregoing, it is assumed that the clock edges corresponding to the two 
signals imder test are separated no more than one clock period from one another. When 
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the corresponding clock edges are separated more than one period, the deterministic 
clock skew value is determined by the difference between the initial phase angles and by 
the sum of the offset times of those clock edges. 



The clock signal distributed from the clock signal source has a close relationship 
with the clock signal source. As a result, the phase noise (timing jitter sequences) in the 
distributed clock signals generally show a similar trend with the phase noise (timing jitter 
sequences) in the clock signal source. Due to this, the timing jitter sequences of the 
plural clock signals distributed from the same clock signal source both show similar 
characteristics to one another (see Figures 5A and 5B). Therefore, the offset amount 
iipffset in the clock edges corresponding to the two signals under test can be determined by 
the correlation between the timing jitter sequences through estimation by finding an 
offset position where the correlation value is the largest. The offset amoxmt noffset in the 
above clock edges can also be determined based on the offset position where the 
correlation value of the instantaneous phase noise becomes the largest 

In addition, the deterministic clock skew can be found by determining the zero- 
crossing timing of each signal under test and calculating the average value of the tune 
difference between the corresponding zero-crossings. 

The clock skew estimation method of the present invention first determines the 
timing jitters A^[n] and A<t)^[n] of the two signals Xj(t), Xk(t) under test as shown in Figure 
4. The determined timing jitter waveforms A<|>^[n] and A(|)^[n] are shown in Figures 5 A 
and 5B respectively. Next, the deterministic skew value x'^^skew between the two signals 
Xj(t), Xk(t) under test is determined. Then, the timmg difference of the clock edges are 
determined by calculating the difference between the timing jitters A^[n] and A^\n] so 
that the random components of the clock skew between signals Xj(t), Xk(t) under test can 
be determined. By calculating the sum of the random components and deterministic 
components in the clock skew, the clock skew T^'^skewW between the signals under test 
can be determined. The determined clock skew T^'^'skewM is shown in Figure 6. 
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Next, the RMS value and the peak-to-peak value in the clock skew will be 
measured from the clock skew sequence T*'^skew[n]. The RMS value T^*''skew,RMs of the 
clock skew is the standard deviation for clock skew T*'\kew[n], and is determined by the 
following equation: 




[sec] (22) 

Here, N is a number of samples of the measured clock skew data. Also, the peak-to-peak 
value T^'^skew,pp of the clock skew is a difference between the maximimi and minimum 
value of T^'^skewM, and is determined by the following equation: 

[sec] (23) 

Figure 7 shows the histogram of the clock skew measured by this clock skew 
measurement method. 

This clock skew estimation method can also measure the clock skews between the 
clock signals with different frequencies. Here, the clock distribution network shown in 
Figure 8 is considered. The clock source PLLg multiphes the system clock CLKq 
provided from an external source by M times and distributes the clocks CLKj and CLKk 
to the network. Figure 9(a) shows the system clock CLKq, and Figure 9(c) shows the 
multiplied clock CLKj. The A0[1] [rad] of the system clock CLKo indicates the timing 
fluctuation between the edge and the ideal clock edge. Therefore, as shown in Figure 
9(b), assuming that the ideal clock edge is multipKed by M times, AGltn/Mj] and A^[n] 
become 1 : 1 when Ae[l] is copied by M-1 times. Here, the LxJ indicates a maximum 
integer where it does not exceed x. By calculating the clock skew between CLKj and 
CLJCg with equation (6), equation (24) is obtauied as follows: 

[sec] (24) 

The deterministic clock skew value T^^skew between the clocks CLKj and CLKg is 
expressed as the time difference between the ideal clock edge (nMT)j of the clock CLKj 
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and the ideal clock edge (nMT)G of the system clock CLKg, and can be determined from 
the initial phase angle of each clock signal as follows: 
T^^={nMT)j-{nMT)^ 
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Here, since the clock CLKj is a clock which is multiplied the system clock CLKq by M 
times, the fundamental period To of the system clock CLKg is equivalent to M times of 
the fundamental period Tj of CLKj (To = MTj). 

In addition, by using an apparatus capable of measuring two channels at the same 
time to first sample CLK^ and CLKg then sample CLKk and CLK^, this clock skew 
measurement method can measure the clock skew between CLKj and CLK^. 

First, the clocks CLKj and CLKg are sampled at the same time to determine the 
skew of the clocks CLKj and CLKg by using equation (6): 



^Skew 



^Skew ^ 
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[sec] (26) 

Next, the clocks CLKk and CLKg are sampled at the same time to determine the skew of 
the clocks CLKk and CLKg in a similar manner: 



^ikew \p\ ~ ^Skew 



[«(^)-A,.["( 



2k 



[sec] (27) 

Lastly, by detennining the difference between the clock skew sequences calculated 
above, the clock skew between CLKj and CLKk can be obtained: 

T^b] = T^[n]-TSi[n] 



-4 
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This clock skew measurement method can not only estimate the clock skews 
between the distributed MPU clock signals mentioned above, but it can also be applied to 
estimate the clock skews of other signals. 

Next, the probability estimating method for peak-to-peak clock skews in the 
present invention will be explained. 

Figure 10 shows the power spectra determined by using the Fast Fourier 
Transform on the clock waveform of the microprocessor. The upper diagram shows the 
quiet mode in the microprocessor, namely the power spectra in an inactive state of the 
microprocessor, and the lower diagram shows the noise mode in the microprocessor, 
namely the power spectra in an active state of the microprocessor. In the quiet mode, 
only the PLL (phase-locked loop) circuit operates to output the clock signals, and is in the 
condition where the clocks are not affected by surrounding circuit operations. In the 
noise mode, all of the L2 (level 2) memories, system bus, core bus, and branch prediction 
units in the microprocessor operate, and is in the condition where the clock is heavily 
affected by the surrounding circuit operations. In either conditions, the line spectrum in 
the clocks is observed at 400MHz where the random instantaneous phase noise is seen in 
a nearby frequency band with respect to the center frequency 400MHz. This indicates 
the existence of narrow-band random data. In addition, the probability density function 
of the timing jitters in the above clock signals is in a Gaussian distribution as shown in 
Figure 11. Therefore, the timing jitter sequences in the clock signals is a Gaussian-based 
random process. 

As mentioned above, the clock skew T^'^skew between the two clock signals is 
expressed as follows. 



and their random component T-'^s is expressed as the difference between the timing jitter 
sequences as follows: 



Skew 
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Therefore, when the probability density function of the timing jitters A^[n] and A^\n] in 
each clock signal shows an average value 0 and a dispersion (variance) a in the Gaussian 
distiibution, 



exp 



I 2^^ J 



(31) 



5 the probability density function of V^^s is expressed as the convolution. 
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where it becomes a Gaussis distribution by central hmit theorem: 



• = V2c 
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(33) 



(34) 



In other words, the random component T'-'^rs in the clock skew is also a Gaussian-based 
random process. 

In the narrow-band random process {Z[n] } , when a certain instantaneous value 
Z[n] is subject to Gaussian distribution, the peak value collection, namely the maximum 
value collection of Z[n] {max(Z[n])}, becomes closer to Rayleigh distribution when the 
degrees of freedom n (number of sample) is made larger. This principle is described, for 
example, page 542 of J. S. Bendat and A. G. Piersol "Random Data: Analysis and 
Measurement Procedure", Edition, published by John Wiley & Sons, Inc. in 1986, or 
pages 90-92 of D. E. Newland ''An Introduction to Random Vibrations, Spectral & 
Wavelet Analysis", published by Longman Scientific & Technical in 1993. 

As explained above, by the existence of the random components of the 
clock skew and that the random components T-^^s being subject to the Gaussian 
distribution, the peak value collection {Zp} = {max (T^'^rsM)} of the random 
components of the clock skew becomes subject to the Rayleigh distribution. The 
Rayleigh distribution is described on pages 30-31 of S. M. Kay "Fundamentals of 
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Statistical Signal Processing: Detection Theory", published by Prentice-Hall, Inc. in 
1998. 

The probability density function Pr(Zp) of the Rayleigh distribution is known to be 
expressed by the following equation: 



f 7^ \ 



exp 



2al 



= 0 



(35) 



Here, Cz is an RMS value of clock skew T''''rs, and is dispersion. The 

Rayleigh probability density fimction is, as shown in Figure 12, Pr(Zp) ^ 0 when Zp > 0. 

Further, when the peak value Zp is subject to the Rayleigh distribution, the 
probability of Zp to become larger than a certain value Zpk is known to be expressed by 
1 0 the following equation (right-tail probability): 



f 72 \ 



(36) 



Also, the standard deviation of Zpk is expressed by the following equation: 



c 7 = . 



4-;r 



(37) 

The probability Pr(Zp > Zpk) is shown in Figure 13. 

Therefore, by setting Zpk as the worst case peak value of the random components 
in the clock skew and measuring the root mean square in the clock skew of the signal 
under test, the probability of the random components in the clock skew of the signal 
under test that exceeds the worst case peak value Zpk can be estimated, where the 
reliability of that clock distribution network becomes higher as the probability becomes 
lower. 

Based on the consideration described above, the probability estimating method for 
peak-to-peak clock skews in the present invention determines the generation probabihty 
of the peak value in the clock skews between the signals under test. 
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Further, by removing the low frequency components from the instantaneous phase 
noise, the probabiUty density function in the timing jitters can become closer to the 
Gaussian distribution, thereby improving the accuracy in the probability estimation. 

Also, when the probability for peak value Zp of the timing jitter of flie input signal 
that exceeds the peak value Zpk is given by the equation (35), the probability for the peak- 
to-peak value Jpp of the jitter that exceeds Zpp as shown below can be obtained from the 
product of the probability for positive peak value to exceed +Zpp/2, and the 
probabiUty for negative peak value Zp" tiiat exceeds -Zpp/2: 
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Based on the foregoing observation, the probabiUty estimation method for peak-to-peak 
clock skews of the present invention determines the generation probability of the peak-to- 
peak value in the clock skews between the signals under test. 

In Figure 14, the samples and peak-to-peak values are plotted for the clock skews 
between the two signals under test, which have been distributed within the 
microprocessor. The upper diagram shows the test results in the quiet mode in the 
microprocessor, and the lower diagram shows the test results in the noise mode. The 
theoretical curve in the diagrams are calculated from the inverse probabiUty Pr(Zpp>Zpp) 
expressed by the equation (38). The test data complies very well with the theoretical 
curves of the Rayleigh distribution (especially in the noise mode). 

Next, the timing jitter estimation method of the present mvention is described. A 
jitter-free clock signal is a square wave with a fundamental frequency f^. These signals 
can be separated into harmonics frequencies including 5), 3fo, 5fo and so forth based on 
Fourier analysis. Since the jitter corresponds to the fluctuation of the fundamental 
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frequency of the signal under test, only the signal components closest to the fundamental 
frequency will be considered m the jitter analysis. 

The fundamental sinusoidal wave componait in the clock signal (signal under 
test) with jitter is e^qjressed as follows, where A represents an ampUtude of the clock 
5 signal and To represents the fundamental period of the clock signal: 

fin ^ 
x(t) = Acosi^{t)) = Acos — t+^o -^{t) 

(39) 

Here, (|)(t) represents an instantaneous phase of the signal under test, and is expressed as a 
sum of the linear phase component 27rt/To having the fundamental period To, an initial 

II phase angle ^o (can be zero in the calculation), and the instantaneous phase noise 

IS: , 

13 10 component A(|)(t). 

p When the instantaneous phase noise component A(|)(t) is zero, an interval between 



pi 



the rising zero-crossing points of the signal under test is merely the constant period To. 
The non-zero instantaneous phase noise component A<|»(t) causes to fluctuate the zero- 
crossing points of the signal under test. In other words, A(|)(nTo) in the zero-crossing 
|^.| indicates the time fluctuation of the zero-crossing points, and is called a timing 

jitter. Therefore, by estimating the instantaneous phase <|>(t) of the signal under test and 
finding the difference between that instantaneous phase and linear phase in the zero- 
crossing points (which corresponds to the phase waveform of the ideal jitter-free clock 
signal), namely the instantaneous phase noise A<|>(t), the timing jitters in the signals under 
20 test can be determined. 

First, the timing jitter estimation method of the present invention converts the 
signal x(t) under test shown in Figure 15 into an analytic signal z(t) of a complex number. 
The converted analytic signal z(t) is shown in Figure 16. In Figure 16, the continuous 
Une indicates the real part of the analytic signal and the broken line indicates the 
25 imaginary part ofthe analytic signal. Next, an instantaneous phase (|)(t) of the signal x(t) 
under test is estimated from the analytic signal z(t). The estimated instantaneous phase 
waveform ^(t) is shown in Figure 17. Then, by conducting a Unear Une fitting for the 
instantaneous phase waveform data based on the minimum square method, and 
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determining the linear phase (|)ii„ear(t) corresponding to the instantaneous phase waveform 
in the ideal jitter-free signal, and calculating the difference between the instantaneous 
phase (|>(t) and that linear phase <|>iinear(t), the instantaneous phase noise A<|>(t) in the signal 
under test is determined. The instantaneous phase noise waveform A<|)(t) determined in 
this manner is shown in Figure 1 8. Further, by sampling the instantaneous phase noise 
waveform A^(t) at a timing closest to each zero-crossing point (analogous zero-crossing 
point) in the real part x(t) of the analytic signal z(t), the instantaneous phase noise in the 
zero-crossing timing nTo, namely the timing jitter A<|>[n] (=A<t>(nTo)), is estimated. The 
estimated timing jitter waveform A<|)[n] is shown in Figure 19. 

With the use of a waveform clipper, the timing jitter estimation method of the 
present invention can estimate the timing jitter with hi^ accuracy by removing the AM 
(amplitude modulation) components from the signals under test and retaining only the 
PM (phase modulation) components corresponding to the jitter. 

Also, by using means for removing low frequency components, the timing jitter 
estimation method of the present invention can remove the low frequency components 
from the instantaneous phase noise signals. 

The method of estimating the instantaneous phase with use of the analytic signal 
is described here. The analytic signal z(t) of the real signal x(t) is defined by a complex 
signal in the following equation: 

z{t)^x{t)+jx(t) 

(40) 

Here, j is an imaginary unit, and an imaginary part x(t) of the complex signal z(t) is a 
Hilbert-transformation of the real part x(t). 

Meanwhile, a ffilbert transformed time wavform x(t) is defined by the following 
equation: 

x(t)=H[x{t)]=^r^dT 

(41) 

Here, x(t) is the convolution of tiie ftmction x(t) and (l/nf). In otiier words, the Hilbert 
tiansformation is equivalent to the x(t) that has passed through a band-pass filter. 
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However, the spectra component in the output x(t) at this time will not change in the 
amplitude, but shift the phase by nil. 

The analytic signal and Hilbert transformation are described, for example, in A. 
Papoulis ''Probability Random Variables and Stochastic Processed', Edition, 
published by McGraw-Hill Book Company, 1984. 

The instantaneous phase waveform <j)(t) of the real signal x(t) is determined by 
using the following equation from the analytic signal z(t): 

At). 

(42) 

Next, the algorithm for estimating the instantaneous phase with use of the Hilbert 
transformation is explained below. First, by applying the Hilbert transfoimation to the 
signal under test shown in Figure 20 



(j)(/)=tan"* 



(O=^cos —t+^^-A^{t)\ 



(43) 



arid determining a signal corresponding to the imaginary part of the complex signal, 



x(t)=H[x{t)]=Asm ^/+(j,„-A(|)(?) 



(44) 



the signal x(t) under test is transformed into an analytic signal in the equation (45): 



z(t) = x(t)+ j'x(t) =Acos — r +<|)o - A<|> (f) 



+ JAsia 
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— f+<|)o-A(j)(f) 



(45) 

The transformed analytic signal is shown in Figure 21. Here, a band-pass filter process is 
^lied to the obtained analytic signal. This is for handling only the signal components 
closest to the fundamental frequency of the signal under test in the jitter analysis because 
the jitter corresponds to the fluctuation of the ftmdamental frequency of the signal under 
test. Next, the phase function (|>(t) is estimated by using equation (42) from the 
determined analytic signal z(t): 
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(46) 



Here, ^(t) is expressed by using the principal value of a phase range between -n and +7r, 
and carries discontinuous points near the conversion point from +n to -ti. The estimated 
phase function <|)(t) is shown in Figure 22. Lastly, by unwrapping the discontinuous 
phase function <j)(t) (namely, appropriately adding the integer multiple of 27i to the 
principal value ^(t)), the discontinuity can be removed to obtain the continuous 
instantaneous phase (j>(t). 



The phase unwr^ping method is described in Donald G. Childers, David P. Skinner, and 
Robert C. Kemerait "JJte C^stnim: A Guide to Processing, volume 65, published by 
Proceedings of IEEE, 1977. The unwrapped continuous instantaneous phase function 
<|)(t) is shown in Figure 23. 

The transformation from the real signal to the analytic signal can be fulfilled by a 
digital signal processing using Fast Fourier Transformation (FFT). 

First, FFT is appUed to the digitized signal x(t) under test as shown in Figure 24 
to obtain both-side power spectra (with positive and negative frequency) X(f) of the 
signal under test. The obtained both-side power spectra X(f) is shown in Figure 25. 
Then, the entire data is emptied except for the data closest to the fundamental frequency 
from the positive frequency components of the spectra X(f), and doubles the positive 
frequency components. Such processes in the frequency domain corresponds to the 
band-limit process of the signal under test in the time domain to obtain the analytic 
signal. The obtained signal Z(f) in the frequency domain is shown in Figure 26. Lastiy, 
by applying an inverse FFT to the obtained signal Z(f), the band-limited analytic signal 
z(t) can be obtained. The band-limited analytic signal z(t) is shown in Figure 27. 

The analytic signal transforaiation using FFT is described, for example, in J. S. 
Bendat and A. G. Piersol "Random Data: Analysis and Measurement Procedure", 
Echtion, pubhshed by John Wiley & Sons, Inc. in 1986, 




^+(j)o-A(|)(r) 
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Alternatively, when the instantaneous phase estimation is the object, the process 
of doubling the positive frequency component can be omitted. 

Next, the detection method for the analogous zero-crossing points will be 
explained. First, a signal value ¥50% at the 50% level of the analytic signal of the input 
signal to be tested is calculated as a zero-crossing level, where the maximum value for 
the real part x(t) of the analytic signal is at a 100% level thereof and the minimum value 
is at a 0% level thereof. Then, the difference between each adjacent sample value and the 
50% level V500/0 , namely, (x(j-l)-V5o%) and (x(j)-V5o%) is calculated to further determine 
the product (x(j-l)-V5oo/o) x (xO)-V5o%). When x(t) is crossing the 50% level, namely the 
zero-crossing level, the symbol of the sampled values (x(j-l)-V5oo/o), (xijyYsm) change 
from negative to positive or vice versa, therefore, when the above product is negative, it 
means that x(t) crosses the zero-crossing level, where a time j-1 or j whichever smaller in 
the absolute value of the sampled values (x(j-l)-V5o%), (xO)-V5oo/o) at that point is 
determined as the analogous zero-crossing point. The waveform of the real part x(t) of 
the analytic signal is shown in Figure 28. The circle marks in Figure 28 show the points 
closest to tiie detected rising zero-crossing point (analogous zero-crossing point). 

Waveform chpping means is able to remove the AM components from the input 
signal to retain only the PM components in the input signal corresponding to the jitters. 
The waveform clipping for the analog or digital input signal is conducted by (1) 
multiplying the signal value by a constant number, (2) replacing the signal value larger 
than the predetermined threshold value 1 with the threshold value 1, and (3) replacing the 
signal value smaller than the predetermined threshold value 2 with the threshold value 2. 
Here, the threshold value 1 is assumed to be larger than the threshold value 2. The clock 
signal with the AM component is shown in Figure 29. From the fluctuation of the time 
waveform envelope in this drawing, the existence of the AM components is obvious. The 
clock signal clipped by the chpping means is shown in Figure 30. Since the time 
waveform shows a constant envelope, the removal of the AM components can be 
confirmed. 

The foregoing is the theoretical aspect of the present invention. Now the 
embodiments of the present invention will be e:q>lained below. 
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Figure 31 shows an example of structure in the probabiUty estunating apparatus 
for peak-to-peak clock skews in the embodiment of the present invention. The 
probability estimating apparatus for peak-to-peak clock skews 100 is comprised of a 
clock skew estimator 101 for estimating the clock skew sequences between a plurality of 
clock signals under test, and a probability estimator 102 for determining and outputting 
the generation probability of the peak-to-peak value in the clock skews between the 
signals imder test based on the clock skew sequences. Further, the probabiHty estimator 
102 is comprised of an RMS (root mean square) detector 103 for detecting the RMS 
value of the jitter waveforms suppUed thereto, a memory 104 for storing a predetermined 
peak-to-peak value, and a probability calculator 105 for calculating the probability of the 
peak-to-peak jitters which exceeds the above predetermined value based on the 
predetermined peak-to-peak value and the detected RMS value. The clock skew 
estimator and its detailed configuration will be explained later. 

Next, with use of the probability estimating apparatus for peak-to-peak clock 
skews 100 in the present invention, the operation to estimate the probability of the peak- 
to-peak clock skews between the plurality of clock signals under test that exceeds the 
predetermined value will be explained. Figure 32 shows the operational procedure of the 
probability estimating method for peak-to-peak clock skews in the present invention. 
First, in step 201, the predetermined peak-to-peak value in the memory 104 is initialized. 
Next, in step 202, the clock skew estimator 101 estimates the clock skew sequences 
among the plurality of clock signals under test and supplies the estimated clock skew 
sequences to the probability estimator 102. Then, in step 203, the RMS detector 103 
calculates the root mean square value in the clock skew sequences supplied from the 
clock skew estimator 101. Lastly, in step 204, the probability calculator 105 calculates 
the probability that the peak-to-peak value in the clock skews between the signals under 
test exceeds the above predetermined value with use of the peak-to-peak value stored in 
the memory 104 and the RMS value, and the process ends. In the step 204 noted above, 
the probability calculator 105 determines the probability of the peak-to-peak value in the 
clock skews exceeding the above predetermined value by using the equation (38), where 
the RMS value obtained from the step 203 is az and the peak-to-peak value stored in the 
memory 104 is Zpp. 
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Alternatively, the set value for determining the probability of the peak-to-peak 
value in the clock skews exceeding the predetermined value can be 2K times (K is 
constant) the RMS value of the clock skews. In this situation, instead of the memory 
104, by providing a constant number multiplier for multiplying the RMS value az 
obtained by the RMS detector 103 by 2K times, the obtained 2Kaz can be input to the 
probability calculator 105 as Zpp, 

The probability estimating apparatus for peak-to-peak clock skews shown in 
Figure 31 can also be structured as an apparatus for determining the probability of the 
peak value of the clock skews exceeding the predetermined value. In this situation, the 
memory 104 stores the predetermined peak value, and the probability calculator 105 
determines the probability of the peak value of the clock skews exceeding the 
predetermined value by using the equation (36). 

Similarly, in the probability estimating method for peak-to-peak clock skews 
shown in Figure 32, by replacing the step 204 for calculating the generation probability 
of the peak-to-peak value in the clock skews between the signals under test with a step of 
calculating the generation probability of the peak value in the clock skews, the peak value 
in the clock skews can be used in the procedure for determining the generation 
probabihty of the peak value in the clock skews exceeding the predetermined value. 

Further, when determining the probability for the peak value in the clock skews 
exceeding the predetermined value, this predetermined value can be K times (K is 
constant) the RMS value in the clock skews. Under this condition, instead of the memory 
104, by providing a constant number multiplier for multiplying the RMS value az 
obtained by the RMS detector 103 by K times, the obtained Kaz can be input to the 
probability calculator means 105 as the peak value Zpk. 

Also, the probability estimating apparatus for peak-to-peak clock skews shown in 
Figure 31 can be structured as an apparatus for determining the generation probability of 
the peak-to-peak value in the clock skews by changing the configuration of the 
probability estimator. Figure 33 shows such an example of structure in the probability 
estimating apparatus for peak-to-peak clock skews used in another embodiment of the 
present invention. This probability estimating apparatus 300 for peak-to-peak clock 
skews in Figure 33 is similar to the probability estimating apparatus for peak-to-peak 
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clock skews shown in Figure 31 except for the probabiHty estimator 102 being replaced 
with a probability estimating apparatus 102a, where it is comprised of an RMS detector 
103 for determining the RMS value of the clock skew sequences supplied thereto, a peak- 
to-peak detector 301 for determining the peak-to-peak value based on the supplied clock 
5 skew sequences, and a probabiUty calculator 105 for calculating the probability of the 
clock skews between the signals under test exceeding the above peak-to-peak value based 
on the peak-to-peak value and the RMS value (for simplicity of explanation, descriptions 
for the duplicated parts are omitted). 

Next, with use of the probability estimating apparatus for peak-to-peak clock 
10 skews 300 in the present invention, the operation for estimating the probability of the 
clock skews between the signals under test exceeding the peak-to-peak value will be 
explained. Figure 34 shows such an operational procedure of the probability estimating 
method for peak-to-peak clock skews in the present invention. This probabiUty 
estimating method for peak-to-peak clock skews is similar to the probability estimating 
15 method for peak-to-peak clock skews shown in Figure 32 except for Ihe removal of the 
rl step 201 for initializing the predetermined peak-to-peak value in the memory, and the 



m 



p f additionof a Step 401 in which, after the RMS vahie in the clock skew sequ^ 



;| between flie signals under test is determined, the above noted peak-to-peak detector 301 

■ calculates the difference between the maximum value and the minimum value in the 

20 clock skew sequences to determine the peak-to-peak value (for simplicity of explanation, 
descriptions for the duplicated parts are omitted). Here, in step 204, the probability 
calcxilator 105 determines the probability of the clock skews exceeding the peak-to-peak 
vialue by using the equation (38), where flie RMS value obtained from the step 203 is az 
and the peak-to-peak value also obtained from the above step 203 is Zpp. Further, Ihe step 
25 203 for determining the RMS value in the above clock skew sequences and the step 401 
for determining the peak-to-peak value are independent from one another, and can 
conduct the procedure in parallel or different order. 

Also, the probability estimating apparatus for peak-to-peak clock skews shown in 
Figure 33 can be structured as an apparatus for determining the generation probabiUty of 
30 the peak value in the clock skews. In this situation, the probability calculator 105 
determines the probability of Ihe peak value in the clock skews exceeding the above 
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noted predetermined value by using the equation (36), and the peak-to-peak detector 301 
determines the maximum or minimum value in the clock skews. 

Similarly, in the probability estimating method for peak-to-peak clock skews 
shown in Figure 34, by replacing the step 204 for calculating the generation probability 
of the peak-to-peak value in the clock skews between the signals under test with a step of 
calculating the generation probability of the peak value in the clock skews, and by 
replacing the step 401 for determining the peak-to-peak value in the clock skews with a 
step for calculating the maximum or minimum value in the clock skew sequences to 
determine the peak value, the procedure can determine the generation probabiUty of the 
peak value in those clock skews. 

Figure 35 shows an example of structure in the clock skew estimator used in the 
embodiment of the present invention. The clock skew estimator 500 is comprised of 
timing jitter estimators 501a and 501b for estimating the timing jitter sequences of the 
clock signals under test, a deterministic clock skew estimator 502 for estimating the 
timing errors between the ideal clock edges of the clock signals under test to determine 
the deterministic components of the clock skews, and a clock skew calculator 503 for 
receiving the above timing jitter sequences to calculate the timing differences and 
outputting the clock skew sequences. In addition to the timing jitter sequences, the 
timing jitter estimators 501a and 501b also estimate the initial phase angles of the signals 
imder test, and output them to the deterministic clock skew estimator 502. The detailed 
structure of the timing jitter estimator will be explained later. 

Next, with the use of the clock skew estimator 500 in the present invention, the 
operation for conducting the clock skew estimation between the clock signals under test 
will be eaplained. Figure 36 shows the procedure for the clock skew estimating method 
in the present invention. First, in step 601, the timing jitter estimators 501a and 501b 
estimate the initial phase angles of the signals under test and the timing jitter sequences. 
Next, in step 602, the deterministic clock skew estimator 502 calculates the difference 
between the initial phase angles of the signals under test received from the timing jitter 
estimators 501a and 501b to estimate the deterministic components of the clock skews 
between the signals under test Lastly, in step 603, the clock skew calculator 503 
estimates the clock skew sequences between the signals under test based on the timing 
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jitter sequences obtained from the timing jitter estimators 501a and 501b and the 
deterministic components of the clock skews obtained from the detemiinistic clock skew 
estimator 502, and the process ends. 

In the step 602, in which the deterministic components of the clock skews 
between the signals under test are estimated, the deterministic clock skew estimator 502 
obtains the deterministic components in the clock skews between the signals und^ test 
by using the equation (16), Also, in Ihe above step 602, the deterministic clock skew 
estimator 502 can determine the absolute value of the equation (16) when necessary. 
Further, in the step 603, where the clock skew sequences between the signals under test 
are estimated, the clock skew calculator 503 determines the clock skew sequences 
between the signals under test by using the equation (6). Moreover, the step 601 , in 
which the initial phase angles of the signals under test and timing jitter sequences are 
estimated, can be replaced with a procedure shown in Figure 44. In addition, the step 
602, in which the deterministic components of the clock skews between the signals under 
test are estimated, can be replaced with a procedure shown in Figure 38. 

The clock skew estimator shown in Figure 35 can be structured as an ^paratus 
for estimating only the random componaits of the clock skews. In this situation, the 
deterministic clock skew estimator 502 for deteraiining the deterministic components of 
the clock skews can be omitted. Similarly, the clock skew estimating method shown in 
Figure 36 can estimate only the random components of the clock skews. In this 
arrangement, the step 602 for estimating the deterministic components of the clock skews 
from the initial phase angles of the signals under test can be omitted. 

The above deterministic clock skew estimator 502 can be implemented by the 
configuration shown in Figure 37. Figure 37 shows an example of structure in the 
deterministic clock skew estimator used in the embodiment of the present invention. The 
deterministic clock skew estimator 700 is comprised of an offset estimator 701 for 
receiving the initial phase angles and the timing jitter sequences of the signals under test 
for estimating the offset between the clock edges corresponding to the signals xmder test 
based on the timing jitter sequences, and a deterministic clock skew calculator 702 for 
calculating the deterministic components of the clock skews between the signals under 
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test based on the initial phase angles and the offset of the clock edges estimated by the 
offset estimator 701 . 

Next, with use of the deterministic clock skew estimator 700 in the present 
invention, the operation for estimatmg the deterministic components of the clock skews 
between the signals under test will be explained. Figure 38 shows the procedure for the 
deterministic clock skew estimating method in the present invention. First, in step 801, 
the offset estimator 701 determines the offset position where the correlation coefficient 
between the timing jitter sequences becomes the largest based on the timing jitter 
sequences of the signals under test, thereby estimating the offset no^et between the 
corresponding clock edges. Then, in step 802, the deterministic clock skew calculator 
702 calculates the deterministic components of the clock skews between the signals under 
test based on the initial phase angles and the offset in the clock edges estimated by the 
offset estimator 701, and the process ends. In the step 802, where the deterministic 
components of the clock skews between the signals under test are calculated, the 
deterministic clock skew calculator 702 obtains the deterministic components of those 
clock skews between the signals under test by using the equation (21). 

Figure 39 shows another structural example of the clock skew estimator used in 
the embodiment of the present invention. The clock skew estimator 900 is comprised of 
timing jitter estimators 501a, 501b, 501c, and 501d for estimating the timing jitter 
sequences A<|)^[n], A(I)^[n], A<^\nl and A(t)^[n] in the clock signals Xj(t), Xg(t), Xk(t), and 
Xg(t) under test respectively, deterministic clock skew estimators 502a and 502b for 
estimating the timing errors between the ideal clock edges of the clock signals Xj(t) and 
Xg(t) under test as well as between the clock signals Xk(t) and Xg(t) under test to estimate 
the deterministic components T^skew and x^^skew of the clock skews, clock skew 
calculators 503a and 503b for receiving the timing jitter sequences and calculating their 
timing differences therebetween to output the clock skew sequences T^**skew[n] and 
T^'^skewW, and another clock skew calculator 901 for receiving the above clock skew 
sequences from the clock skew calculators 503a and 503b to determine the difference 
between the clock skew sequences and estimate the clock skew sequence T^'^skewM- For 
simplicity of explanation, description for the duplicated parts is omitted. 
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Next, with use of the clock skew estimator 900 in the present invention, the 
operation for conducting the clock skew estimation among the signals under test will be 
cjqplained. Figure 40 shows the procedure of the clock skew estimation method in the 
present inveution. First, in step 1001, timing jitter estimators 501a and 501b estimate the 
initial phase angles <|»'o and and timing jitter sequences A<t>' [n] and A^^[n] of the 
signals Xj(t), Xg(t) under test. Then, in step 1002, the deterministic clock skew estimator 
502a calculates the difference between the initial phase angles and of the signals 
under test obtained from the timing jitter estimators 501a and 501b to estimate the 
deterministic component 'r^''skew of the clock skews between the signals under test. Then, 
in step 1003, the clock skew calculator 503a estimates the clock skew sequence T^''skew[n] 
between the signals under test based on the timing jitter sequences A<}x'[n] and A^\n] 
obtained from the above timing jitter estimators 501a and 501b and based on the 
deterministic component T^''skew of the clock skews obtained from the deterministic clock 
skew estimator 502a. 

Next, in step 1 004, the timing jitter estimators 501 c and 501 d estimate the initial 
phase angles <|)% and <t>^o and the timing jitter sequences A<^\n] and A{|)^[n] of the signals 
Xk(t), Xg(t) under test. In step 1005, the deterministic clock skew estimator 502b 
calculates the difference between the initial phase angles ^\ and of tiie signals under 
test obtained from the timing jitter estimators 501c and 501d to estimate the deterministic 
component T^stew of the clock skews between the signals under test. Then, in step 1006, 
the clock skew calculator 503b estimates the clock skew sequence T^^skewM between tiie 
signals under test based on tiie timing jitter sequences A<^\n] and A<|)^[n] obtained from 
the above timing jitter estimators 501c and 501d and based on the deterministic 
component x^^'skew of the clock skews obtained from the deterministic clock skew 
estimator 502b. 

Lastly, in step 1007, the clock skew calculator 901 estimates tiie clock skew 
sequence T-i^skewEn] between flie signals Xj(t), Xk(t) under test based on the clock skew 
sequences T^skew[n] and T^-^skewM obtained from tiie clock skew calculators 503a and 
503b, and tiie process ends. In tiie above step 1007, in which tiie clock skew sequences 
between tiie signals Xj(t), Xk(t) under test are estimated, tiie clock skew calculator 901 
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determines the clock skew sequences between the signals under test by using the equation 
(28). To simplify the explanation, descriptions for the duplicated parts are omitted. 

The clock skew estimator shown in Figure 39 can be structured as an apparatus 
for estimating only the random components of the clock skew. Under such a 
5 construction, the deterministic clock skew estimators 502a and 502b for determining the 
deterministic components of the clock skews can be eliminated. Similarly, the clock 
skew estimating method shown in Figure 40 can estimate only the random components of 
the clock skew. Under such a procedure, the steps 1002 and 1005 for estimating the 
deterministic components of the clock skews from the initial phase angles of the signals 
10 under test can be eUminated. 

Figure 41 shows another structural example of the clock skew estimator used in 
f I the embodiment of the present invention. This clock skew estimator 11 00 is similar to 

fU clock skew estimator shown in Figure 35 except for the inclusion of a frequency 
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multiplier 1 101 for receiving the timing jitter sequences estimated by the timing jitter 
U 15 estimator 501b to copy the timing jitter sequences by, for example, M-1 times, and to 
£; | determine the timing jitter sequences when the signals under test are multiplied by M 

^ times (in order to simplify the explanation, descriptions for the duplicated parts are 

tij omitted). 

Next, with the use of the clock skew estimator 1 100 in the present invention, the 
20 operation for conducting the clock skew estimation in the signals under test will be 
explained. Figure 42 shows another procedure of the clock skew estimation method in 
the present invention. This clock skew estimation method is similar to the clock skew 
estimation method shown in Figure 36 except for having a step 1201 in which, after 
determining the timing jitter sequences, the frequency multiplier 1101 copies the timing 
25 jitter sequences estimated by the timing jitter estimator 5 10b by, for example, M-1 times, 
to determine the timing jitter sequences when the signals under test are multiplied by M 
times (in order to simplify the explanation, descriptions for duplicated parts are omitted). 
Under this situation, in the above step 602 for estimating the deterministic components of 
the clock skews between the signals under test, the deterministic clock skew estimator 
30 502 obtains the deterministic components of the clock skews between the signals under 
tested by using the equation (25). Also, in the above noted step 603 in which the clock 
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skew sequences between the signals under test are estimated, the clock skew calculator 
503 determines the clock skew sequences between the signals under test by using the 
equation (24). 

The clock skew estimator shown in Figure 41 can be structured as an ^aratus 
for estimating only the random components of the clock skews. In such a situation, the 
deterministic clock skew estimator 502, which determines the deterministic components 
ofthe clock skews, can be eliminated. Similarly, the clock skew estimating method 
shown in Figure 42 can also estimate only the random components ofthe clock skews. 
At this time, the step 602 in which the deterministic components ofthe clock skews can 
be estimated from the initial phase angles of the signals under test, can be omitted. 

Further, the above frequency multiplier can be incorporated into the clock skew 
estimator shown Figure 39. In such a construction, the frequency multiplier is inserted 
directly at the output of the timing jitter estimators 501b and 501d. Similarly, the step for 
multiplying the frequency can be added to the clock skew estimating method shown in 
Figure 40. At this time, the step to multiply the frequency is inserted after the steps 1001 
and 1004 in which the timing jitters are estimated. 

Figure 43 shows a structural example ofthe timing jitter estimator used in the 
clock skew estimator of the present invention. The timing jitter estimator 1300 is 
comprised of an analytic signal transformer 1301 for transforming the signal under test 
into a band-limited analytic signal of complex number, an instantaneous phase estimator 
1302 for determining the instantaneous phase of the analytic signal transformed by the 
analytic signal transformer 1301, a linear trend remover 1303 for removing the linear 
phase from the mstantaneous phase estimated by the instantaneous phase estimator 1302 
to obtain the instantaneous phase noise, and a zero-crossing resamplerl304 for receiving 
the instantaneous phase noise estimated by the linear trend remover 1303 to resample 
only the instantaneous phase noise data closest to the zero-crossing timings ofthe real 
piart ofthe analytic signal and to output the timing jitter sequences. The analytic signal 
tiransformer 1301 can use the stiucture shown in Figure 45, 47, or 49. Also, the analytic 
signal fransformer 1301 can be structured to freely change the pass band ofthe signals. 
Further, the linear trend remover 1303 determines the initial phase angles ofthe signals 
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under test as well as the instantaneous phase noise, and outputs them to the deterministic 
clock skew estimator. 

Next, with use of the timing jitter estunator 1300 in the present invention, the 
operation in the method for estimating the initial phase angles as well as the timing jitter 
sequences in the signals under test will be explained. Figure 44 shows the procedure of 
the timing jitter estimation method in the present invention. First, in step 1401, the 
analytic signal transformer 1301 transforms the signals under test into an analytic signal, 
where the predetermined frequency components are selectively passed there through. 
Next, in step 1402, the instantaneous phase estimator 1302 estimates the instantaneous 
phase of the signal under test by using the analytic signal obtained from the analytic 
signal transformer 1301. Then, in step 1403, the linear trend remover 1303 estimates the 
linear phase corresponding to an ideal clock signal based on the instantaneous phase 
estimated by the instantaneous phase estimator 1302, and determines the initial phase 
angle of the signal under test. 

Then, in step 1404, the linear trend remover 1303 removes the Unear phase from 
the instantaneous phase to estimate the instantaneous phase noise. Lastly, in step 1405, 
the zero-crossmg resampler 1304 resamples only the instantaneous phase noise closest to 
the zero-crossing timing of the real part of the analytic signal based on the instantaneous 
phase noise estimated by the linear trend remover 1303 to estimate the timing jitter 
sequences, and the process ends. The above noted step 1401 for transforming the signal 
under test into the analytic signal can be conducted by the procedure shown in Figure 46, 
48, or 50. 

Figure 45 shows a structural example of the analytic signal transformer used in 
the timing jitter estimator 1300 of the present invention. The analytic signal transformer 
1500 is comprised of a band-pass filter 1501 for extracting only the components closest 
to the fundamental frequency of the signal under test to band-hmit the signal under test, 
and a Hilbert transformer 1 502 for Hilbert-transforming the output signal from the band- 
pass filter 1501 to generate a Hilbert transformed pair. The band-pass filter 1501 can be 
an analog filter or a digital filter, or can be implemented through digital signal processing 
such as FFT. Further, the band-pass filter 1501 can be structured to freely change the 
pass band of the signal. 
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Next, with the use of the analytic signal transformer 1500 in the present invention, 
the operation for transforming the signal under test into band-limited analytic signal will 
be explained. Figure 46 shows the procedure for the signal transformation method of the 
present invention. First, in step 1601, the band-pass filter 1501 removes only the 
components closest to the fundamental firequency fi-om the signal under test to band-limit 
those signal under test. Then, in step 1602, the Hilbert transformer 1502 applies the 
Hilbert transformation to the band-limited signal under test, and generates a Hilbert 
transformed pair of the input signal corresponding to the imaginary part of the analytic 
signal. Lastly, in step 1603, the analytic signal transformer 1500 generates the output 
signal from the band-pass filter 1501 as a real part of the analytic signal and the output 
signal fi-om the Hilbert transformer 1502 as an miaginary part of the analytic signal, then 
the process ends. 

Figure 47 shows another structural example of the analytic signal transformer 
used in the timing jitter estimator 1300 of the present invention. The analytic signal 
transformer 1700 is comprised of a time domain to firequency domain transformer 1701 
for transforming the signal under test into both-side spectra signal in a fi-equency domain, 
a bandwidth limiter 1702 for extracting only the components closest to the positive 
fimdamental firequency from the both-side spectra signal in the frequency domain, and a 
frequency domain to time domain transformer 1703 for transforming the output of the 
bandwidth limiter 1702 back into the time domain signal The time domain to frequency 
domain transformer 1701 and the frequency domain to time domain transformer 1703 can 
be implemented by FFT and inverse FFT, respectively. Also, the bandwidth limiter 1702 
can be structured to freely change the pass band of the signal. 

Next, with use of the analytic signal transformer 1700 in the present invention, the 
operation for transforming the signals under test into band-limited analytic signals will be 
explained. Figure 48 shows another procedure of the signal transformation method in the 
present invention. First, in step 1801, the time domain to frequency domain transformer 
1701 applies the FFT to the signal under test to transform the time domain signal into the 
both-side spectra signal in the frequency domain. Then, in step 1802, the bandwidth 
hmiter 1702 replaces the negative frequency components with zero for the transformed 
both-side spectra signal in the frequency domain. Then, in step 1803, the bandwidth 
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limiter 1702 retains only the components closest to the fundamental frequency of the 
signal under test in the one-side spectra signal, where the negative frequency components 
are replaced with zero, to band-limit the signal in the frequency domain. Lastly, in step 
1804, the frequency domain to time domain transformer 1 703 applies the inverse FFT to 
5 the band-limited one-side spectra signal to transform the frequency domain signal into 
analytic signal in the time domain, and then the process ends. The order of the above 
steps 1802 and 1803 can be switched to one another, i.e., after the step of retaining only 
the components closest to the ftindamental frequency of the signal under test to replace 
the remaining frequency components by zero and band-limit signal in the frequency 
10 domain, the step for replacing the negative frequency components in the both-side spectra 
signal with zero can be performed. 

Figure 49 shows a further example of structure in the analytic signal transformer 

IJ used in the timing jitter estimator 1300 of the present invention. The analytic signal 

ik 

|«l . transformer 1 900 is comprised of a buffer memory 1 901 for storing the signal under test, 

15 a signal selector (waveform data selector) 1902 for sequentially extracting the signal from 
F| the buffer memory 1901 while overlapping the part of the exfracted signal with the one 

previously extracted, a window ftmction multiplier 1903 for multiplying a window 
fcl ftmction by each extracted part of the signal, a time domain to fi:^ency domain 

I'll transformer 1 904 for transforming each part of the signal multiplied by the window 

20 ftmction into both-side spectra signal in the frequency domain, a bandwidth limiter 1905 
for extracting only the components closest to the positive fimdamental frequency of the 
signal under test from the both-side spectra signal transformed into the frequency 
domain, a frequency domain to time domain transformer 1906 for transforming the 
output of the bandwidth limiter 1 905 back into the time domain signal, and an inverse 
25 window ftmction multiplier 1907 for multiplying the inverse window ftmction by the 
signal transformed into the time domain to obtain the band-limited analytic signal. The 
time domain to frequency domain transformer 1904 and the frequency domain to time 
domain transformer 1906 can be implemented by using the FFT and the inverse FFT, 
respectively. Further, the bandwidth Umiter 1905 can be structured to freely change the 
30 pass band of the signal. 
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Next, with use of the analytic signal transformer 1900 in the present invention, the 
operation for transforming the signal under test into the band-limited analytic signal will 
be Gxplmnod, Figure 50 shows a further procedure of the signal transformation method 
of the present invention. First, in step 2001, the signal under test is stored in the buffer 
memory 1901 . Then, in step 2002, the signal selector 1902 extracts a part of the signal 
stored in the buffer memory 1901 . Thm, in step 2003, the window function multiplier 
1903 multipUes the window function to the extracted part of the signal. Next, in step 
2004, the time domain to frequency domain transformer 1904 appUes the FFT process to 
the signal multipUed by the window function to transform the time domain signal into the 
botii-side spectra signal in the frequency domain. Then, in step 2005, the bandwidtii 
Umiter 1905 replaces the negative frequency components with zero for the both-side 
spectra signal in the frequency domain. 

Next, in step 2006, the bandwidth limiter 1905 retains only the components 
closest to the fundamental frequency of the signal under test in the one-side spectra 
signal, where the negative frequency components are replaced with zero, to replace the 
remaining frequency components with zero, thereby band-limiting the frequency domain 
signal. Then, in stq) 2007, the frequency domain to time domain transformer 1906 
applies the inverse FFT to the band-limited one-side spectra signal in the frequency 
domain to transform the frequency domain signal into the time domain signal. Then, in 
step 2008, the inverse window function multipKer 1907 multiplies the inverse window 
function produced in the step 2003 by the inverse transformed time domain signal to 
determine the band-limited analytic signals. 

Lastly, in step 2009, inspections are performed to see if there is any unprocessed 
data existing in the buffer memory. If it does exist, in step 2010, the signal selector 1902 
extracts the signal from the buffer memory while overlapping the part of the signal with 
the previously extracted signal, and the steps 2003-2009 are repeated thereafter. If it does 
not exist, the process will end. The order of the above steps 2005 and 2006 can be 
switched to one another, i.e., after the step of retaining only the components closest to the 
fundamental frequency of the signal vmder test by replacing the remaining frequency 
components with zero and band-limit those frequency domain signal, the step for 
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replacing the negative frequency components with the both-side spectra signal with zero 
can be performed. 

Figure 5 1 shows a further structural example of the clock skew estimator used in 
the embodiment of the present invention. This clock skew estimator 2100 is similar to 
5 the clock skew estimator shown in Figure 35 except for the inclusion of AD converters 
12101a and 2101b for digitizing the analog signals under test to be converted into digital 
signals. In order to simplify the explanation, descriptions for duplicated parts are 
omitted. For the AD converter, a high speed AD transformer, digitizer, or digital 
sampling oscilloscope can be preferably used 
10 Next, with use ofthe clock skew estimator 2100 in the present invention, the 

operation ofthe clock skew measurement for the signals under test will be explained* 
;| Figure 52 shows such a procedure of the clock skew estimation method in the present 

1^1 invention. This clock skew estimating method is similar to the clock skew estimating 

^1 method shown in Figure 36 except for the inclusion of a step 2201 in which, prior to the 

|«J 15 Start of the procedure, the AD converters 2101a and 2101b sample (digitize) the analog 
signals under test, which are subject to the clock skew measurement, and convert them 
M ■ into the digital signals. In order to sinipUfy the explanation, descriptions the 
|i J duplicated parts are omitted. 

;|- 1 The above noted AD converters can be incorporated in the clock skew estimator 

20 1 100 having the frequency multiplier shown in Figure 41. In such a situation, the 

procedure ofthe clock skew estimation method shown in Figure 42 additionally includes 
a step where, prior to the start of the procedure, the AD converters sample (digitize) the 
analog signals under test, which are subject to the clock skew measurement, and convert 
them into the digital signals. 
25 Furflier, the above AD transformers can be incorporated in the clock skew 

estimator 900 shown in Figure 39. Similarly, the step for sampling (digitizing) the 
signals under test to transform them into the digital signals can be incorporated into the 
clock skew estimating method shown in Figure 40. 

Figure 53 shows a further structural example of the clock skew estimator used in 
30 the embodiment ofthe present invention. This clock skew estimator 2300 is similar to 
the clock skew estimator shown in Figure 35 except for the inclusion of waveform 
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clippers (waveform clipping means) 2301a and 2301b for removing the AM components 
in the signals under test. For simpUcity of explanation, descriptions for the dupUcated 
parts are omitted. 

Next, with use of the clock skew estimator 2300 in the present invention, the 
5 operation for the clock skew measurement of the signals under test will be explained. 
Figure 54 shows another procedure of the clock skew estimating method of the present 
invention. This clock skew estimating method is similar to the jitter measurement 
method shown in Figure 36 except for the inclusion of a step 2401, where the waveform 
cUppers 2301a and 2301b remove the AM components from tihie signals under test prior 
10 to the start of the procedure. In order to shorten the explanation, descriptions for the 
duplicated parts are omitted. 

The above waveform clippers can be incorporated in the clock skew estimator 
1 100 having the frequency multiplier shown in Figure 41. In such an arrangement, the 
procedure for the clock skew estimation method shown in Figure 42 is provided with a 
15 step where the waveform clippers remove the AM components from the signals under 
test prior to the start of the procedure. 

Also, the above waveform clippers can be incorporated into the clock skew 
estimator 900 shown in Figure 39. Similarly, the step where the waveform clipper 
ft J removes the AM components of the signals under test can be incorporated in the clock 

20 skew estimating method shown in Figure 40. 

Figure 55 shows a ftirther structural example of the timing jitter estimator used in 
the embodiment of the present invention. This timingjitter estimator 2500 is similar to 
the timing jitter estimator shown in Figure 43 except for the inclusion of a low frequency 
component remover 2501 which receives tiie instantaneous phase noise and removes the 
25 low frequency components from the instantaneous phase noise and provides the 
instantaneous phase noise to the zero-crossing resampler. In order to simplify the 
explanation, descriptions for the duplicated parts are omitted. 

Next, with use of the timing jitter estimator 2500 in the present invention, the 
operation for estimating the initial phase angles of the signals under test and timing jitter 
30 sequences will be explained. Figure 56 shows such a procedure for the timing jitter 
estimation method of the present invention. This timingjitter estimating method is 
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similar to the timing jitter estimating method shown in Figure 44 except for the inclusion 
of a step 2601 in which, after estimating the instantaneous phase noise, the low frequency 
component remover 2501 removes the low frequency components from the instantaneous 
phase noise. For simplicity of e}q)lanation, descriptions for the duplicated parts are 
omitted. 

According to the probability estimation apparatus for peak-to-peak clock skews as 
well as the probability estimation method for peak-to-peak clock skews of the present 
invention, by assuming the linear clock skews to be Gaussian-based random processes 
and determining the Rayleigh probability density distribution fimction of the peak-to- 
peak value in the clock skews, the estimation of the probability of the peak-to-peak value 
in the clock skews exceeding the predetermined value, for which no traditionally 
effective methods existed, can be achieved, resulting in the dramatic improvement in the 
effectiveness of the product reliability analysis. 

Furthermore, according to the probability estimation apparatus for peak-to-peak 
clock skews as well as the probabiHty estimation method for peak-to-peak clock skews of 
the present invention, by measuring the peak-to-peak value in the clock skews and 
calculating the generation probability thereof, the generation probability of the peak-to- 
peak value can be examined as to whether it satisfies the product specifications or not, 
thereby dramatically improving the effectiveness of the reliability analysis of products 
which was not possible in the conventional technology. 
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